The hyperfine energy levels of alkali metal dimers: 
ground-state homonuclear molecules in magnetic fields 
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We investigate the hyperfine energy levels and Zeeman splittings for homonuclear alkali-metal 
dimers in low-lying rotational and vibrational states, which are important for experiments designed 
to produce quantum gases of deeply bound molecules. We carry out density-functional theory 
(DFT) calculations of the nuclear hyperfine coupling constants. For nonrotating states, the zero- 
field splittings are determined almost entirely by the scalar nuclear spin-spin coupling constant. By 
contrast with the heteronuclear case, the total nuclear spin remains a good quantum number in a 
magnetic field. We also investigate levels with rotational quantum number TV = 1, which have long- 
range anisotropic quadrupole-quadrupole interactions and may be collisionally stable. For these 
states the splitting is dominated by nuclear quadrupole coupling for most of the alkali-metal dimers 
and the Zeeman splittings are considerably more complicated. 
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I. INTRODUCTION 

There have been enormous advances over the last year 
in experimental methods to produce ultracold molecules 
in their rovibrational ground state at microkelvin tem- 
peratures. Ospelkaus et al. [l| produced KRb mole- 
cules in high-lying states by magnetoassociation (Fes- 
hbach resonance tuning) and then transferred them by 
stimulated Raman adiabatic passage to levels of the 1 E + 
ground state bound by more than 10 GHz. This was 
then extended by Ni et al. 2] to produce molecules in 
(v, N) = (0, 0), where v and N are the quantum numbers 
for molecular vibration and mechanical rotation. Danzl 
et al. 0, 0] have carried out analogous experiments on 
Cs dimers, while Lang et al. [!, [|| have produced Rb2 
molecules in the lowest rovibrational level of the lowest 
triplet state. There have also been considerable successes 
in direct photoassociation to produce low-lying states 

A major goal of the experimental work is to produce 
a stable molecular quantum gas. However, such a gas 
can form only if (i) a large number of molecules are in 
the same hyperfine state and (ii) the molecules are stable 
to collisions that occur in the gas. In particular, inelas- 
tic collisions that transfer internal energy into relative 
translational energy cause heating and/or trap loss. It is 
thus very important to understand the hyperfine struc- 
ture of the low-lying levels and its dependence on ap- 
plied electric and magnetic fields. In a previous paper, 
we explored the hyperfine levels of heteronuclear alkali 
metal dimers in rotationless levels [ll|. The purpose of 
the present paper is to extend this work to homonuclear 
molecules, which have important special features. We 
also explore N = 1 levels, which may be collisionally sta- 



ble for homonuclear molecules and which interact with 
longer-range forces than N = levels. 

II. MOLECULAR HAMILTONIAN 

The Hamiltonian of a diatomic molecule in the pres- 
ence of an external magnetic field can be decomposed into 
five different contributions: the electronic, vibrational, 
rotational, hyperfine and Zeeman terms. For E mole- 
cules in a fixed vibrational level, the first two terms take 
a constant value and the rotational, hyperfine and Zee- 
man parts of the Hamiltonian may be written [ll| [H, [H| 



H — H r , 



H 



hi 



Hz, 



(1) 



where 



H ro t 

Hhf 



B V N 2 — D V N 2 ■ TV 2 ; (2) 

2 



H z = 



Hq + i?iN + H t + H sc — Vi : Qi 

i=l 

2 

* N ■ It + cs h ■ T ■ I 2 + c 4 h ■ h, (3) 

i=l 

2 

-g r [i N N ■ B -J29iVNli ■ B{1- <7i). (4) 



'Electronic address: 
' Electronic address: 



E-mail: Jesus.Aldcgundc@durham.ac.uk 



E-mail: J.M.Hutson@durham.ac.uk 



where the index i refers to each of the nuclei in the 
molecule. N, Jj and I2 represent the operators for me- 
chanical rotation and for the spins of nuclei 1 and 2. The 
rotational and centrifugal constants of the molecule are 
given by B v and D v (but centrifugal distortion is ne- 
glected in the present work). We use N rather than J 
for mechanical rotation because we wish to reserve J for 
the angular momentum including electron spin for future 
work on triplet states. 

The hyperfine Hamiltonian of equation [3] consists 
of four different contributions. The first is the elec- 
tric quadrupole interaction Hq, with coupling constants 
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(eqQ)i and {eqQ)^- It represents the interaction of the 
nuclear quadrupoles (eQi) with the electric field gra- 
dients qi created by the electrons at the nuclear posi- 
tions. The second is the spin- rotation term Hin, which 
describes the interaction of the nuclear magnetic mo- 
ments with the magnetic moment created by the rota- 
tion of the molecule. Its coupling constants are c\ and 
c 2 . For a homonuclear molecule with identical nuclei, 
(eqQ)i — (eqQ) 2 and c\ — c 2 . The last two terms repre- 
sent the interaction between the two nuclear spins; there 
is both a tensor component H t , with coupling constant 
C3, and a scalar component H sc , with coupling constant 
C4. The second-rank tensor T represents the angular part 
of a dipole-dipole interaction. 

The Zeeman Hamiltonian Hz has both rotational and 
nuclear Zeeman contributions characterized by g-factors 
g r , gi and <? 2 . For homonuclear molecules g\ = g 2 . The 
nuclear shielding tensor cr, is approximated here by its 
isotropic part af, terms involving the anisotropy of er.j are 
extremely small for the states considered here. 

The nuclear ^-factors and the qua drup ole moments of 
the nuclei are experimentally known [151 ]. 

For homonuclear molecules we neglect the effect of 
electric fields, though in principle there are small ef- 
fects due to anisotropic polarizabilitics and the molecular 
quadrupole moments can interact with the gradients of 
inhomogeneous fields. 



III. EVALUATION OF THE COUPLING 
CONSTANTS 

The rotational g-factors are known experimentally for 
all the homonuclear alkali metal dimers [l6| . However, 
the only such species for which the nuclear hyperfine cou- 
pling constants have been determined accurately is Na 2 
[l7| . We have therefore evaluated the remaining coupling 
constants using density-functional theory (DFT) calcula- 
tions performed with the Amsterdam density functional 
(ADF) package [HI, [l!| with all-electron basis sets and 
including relativistic corrections. A full description of 
the basis sets, functionals, etc. used in the calculations 
has been given in our previous paper on heteronuclear 
systems [ll|. In the present work, the calculations were 
carried out at the equilibrium geometries (R e = 2.67 A 
for Li 2 [113.08 A for Na 2 (2l|, 3.92 A for K 2 H, 4.21 
A for Rb 2 [Si] and 4.65 A for Cs 2 0). This give results 
that are approximately valid not only for v — states 
but also for other low-lying vibrational states. 

The values for the coupling constants are given in ta- 
ble HI It may be seen that the DFT results for Na 2 are 
within about 30% of the experimental values, and similar 
accuracy was obtained for other test cases in our previ- 
ous work [llj]. The accuracy is likely to be comparable 
for the other cases studied here. This level of accuracy 
is adequate for the purpose of the present paper, which 
aims to explore the qualitative nature of the Zeeman pat- 
terns. Most of our conclusions are insensitive to the exact 



magnitudes of the coupling constants. 

IV. HYPERFINE ENERGY LEVELS 

Our previous work [ll[ showed that the zero-field split- 
ting for heteronuclear diatomic molecules in N = states 
is determined almost entirely by the scalar nuclear spin- 
spin interaction. This remains true for homonuclear mol- 
ecules in TV = states. We show below that for N > 
the electric quadrupole interaction is dominant for all the 
homonuclear dimers except Cs 2 and 6 Li 2 , with smaller 
but significant contributions from the remaining coupling 
constants. 

For all systems except Li 2 , the scalar spin-spin cou- 
pling is considerably stronger than the spin-rotation and 
tensor spin-spin couplings. Knowledge of the nuclear 
quadrupole coupling constant eQq and the scalar spin- 
spin coupling constant C4 is therefore sufficient to un- 
derstand the hyperfine splitting patterns. We will focus 
here on 85 Rb 2 and 87 Rb 2 , which form a convenient pair 
that approximately cover the range of values of the ratio 
\cn/(eqQ)\. Lang et al. [y] have produced 87 Rb 2 in the 
lowest rovibrational level of the lowest triplet state, but 
as far as we are aware not yet in the singlet state. 

The hyperfine energy levels are obtained by diagonal- 
izing the matrix representation of the Hamiltonian (fT]) 
in a basis set of angular momentum functions. In order 
to facilitate the assignment of quantum numbers to the 
energy levels, two different basis sets are employed, 

\{Iih)IMiNM N ) (spin-coupled basis); (5) 
\{hh)INFM F ) (fully coupled basis). (6) 

where / and F are the total nuclear spin and total angular 
momentum quantum numbers and Mj and Mp are their 
projections onto the Z axis defined by the external field. 
The matrix elements of the different terms in the Hamil- 
tonian in each of the basis sets are calculated through 
standard angular momentum techniques 25]. Explicit 
expressions are given in the Appendix. 

For homonuclear molecules, nuclear exchange symme- 
try dictates that not all possible values of the total nu- 
clear spin I can exist for each rotational level. For mol- 
ecules in 1 E + states, only even / values can exist for 
even N and only odd / for odd N. This is true for ei- 
ther fermionic or bosonic nuclei but is reversed for 3 S + 
states. Table [XT] summarizes the I-N pairs compatible 
with the antisymmetry of the wave function under nu- 
clear exchange for the Rb 2 isotopomers. 

A. Zeeman splitting for TV = homonuclear alkali 
dimers 

The Zeeman splittings for the TV = hyperfine levels 
of 85 Rb 2 and 87 Rb 2 are shown in figure [TJ The zero-field 
splittings are in most respects similar to those found for 
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TABLE I: Nuclear quadrupole moment (Q), electric quadrupole coupling constant (eqQ), nuclear g-factor (g), spin-rotation 
coupling constant (ci), tensor spin-spin coupling constant (03), scalar spin-spin coupling constant (04), absolute value of the 
d/(eqQ) ratio, isotropic part of the nuclear shielding (a) and rotational g-factor (g T ) for the homonuclear alkali dimers. All 
the quantities except the nuclear quadrupole moments, the nuclear ^-factors and the rotational (/-factors were evaluated using 
DFT calculations (see section iHl)) . Both experimental [T?} and theoretical results are presented for Na2. 





Q(fm 2 ) eqQ (MHz) 


9 


Cl (Hz) c 3 (Hz) 


c^xiz ) 


\C4/{eqQ)\ 


a (ppm) 




6 Li 2 


-0.082 


0.00123 


0.822 


161 


137 


32 


0.026 


102 


0.1259 


7 Li 2 


-4.06 


0.0608 


2.171 


365 


955 


226 


0.0037 


102 


0.1080 


23 Na 2 (Exp.) 


10.45 


-0.459 


1.479 


243 


303 


1067 


0.0023 




0.0386 


23 Na 2 (DFT) 




-0.456 




299 


298 


1358 


0.0030 


613 




39 K 2 


5.85 


-0.290 


0.261 


35 


5 


106 


0.00036 


1313 


0.0212 


40 K 2 


-7.3 


0.362 


-0.324 


-42 


8 


163 


0.00045 


1313 


0.0207 


41 K 2 


7.11 


-0.353 


0.143 


18 


2 


32 


0.000091 


1313 


0.0202 


85 Rb 2 


27.7 


-2.457 


0.541 


63 


30 


2177 


0.00089 


3489 


0.0095 


87 Rb 2 


13.4 


-1.188 


1.834 


209 


346 


25021 


0.021 


3489 


0.0093 


133 Cs 2 


-0.355 


0.0486 


0.738 


96 


119 


12993 


0.27 


6461 


0.0054 



TABLE II: Values of / permitted by the nuclear exchange 
symmetry for even and odd rotational levels of 85 Rb 2 and 
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I 






even 


0, 2, 4 


85 Rb 2 (/ Rb 


= 5/2) 


odd 
even 


1, 3, 5 
0, 2 


87 Rb 2 (/ Rb 


= 3/2) 


odd 


1, 3 



heteronuclear molecules in the ground rotational state 
[Til ] . The similarities can be summarized as follows: 

• The scalar nuclear spin-spin interaction and the nu- 
clear Zeeman effect are the only two terms in the 
molecular Hamiltonian with nonzero diagonal ele- 
ments for N = 0. 

• The electric quadrupole and the tensor nuclear 
spin-spin interactions are not diagonal in N, cou- 
pling the N, N + 2 and N — 2 rotational levels. 
This means that the energy levels should be con- 
verged by including in the calculations as many ro- 
tational levels as necessary. However, the coupling 
constants eQq and C3 are very much smaller than 
the rotational spacings, so that in practice it is ad- 
equate to include one excited rotational level. Con- 
vergence for N — is reached with N max = 2 and 
convergence for N = 1 is reached with N max = 3. 

• The scalar spin-spin interaction is diagonal in both 
the spin-coupled and fully coupled basis sets, which 



for N = are identical, 

({N = G)IMi\c A h ■ I 2 \(N = Q)IMi) 

= ~c 4 [/(/+l)-2J Rb (J Rb + l)]. (7) 

Except for a very small contribution coming from 
the coupling with N = 2 levels, these diagonal ele- 
ments determine the zero-field splitting. 

Despite the similarity of the zero-field levels, there are 
important differences between the Zeeman splittings for 
heteronuclear and homonuclear molecules. For heteronu- 
clear dimers [ll[ , levels with the same Mi but different I 
exhibit avoided crossings as a function of magnetic field. 
Because of this, I is no longer a good quantum number 
at high field but the individual nuclear spin projections 
Mn and M/ 2 become nearly conserved. For homonuclear 
dimers, however, different energy levels that correspond 
to the same value of Mi are parallel, so that no avoided 
crossings appear as a function of the field. Both I and Mi 
remain good quantum numbers regardless of the value of 
the magnetic field but M/i and M/ 2 are not individually 
conserved. This is illustrated in figure[T] It arises because 
the nuclear Zeeman term, which is the only nondiagonal 
term for N = in the heteronuclear case, is diagonal for 
homonuclear molecules. Its nonzero elements are given 
by equation IA.11I of the Appendix, 

((N = Q)IMi\H IZ \(N = Q)IMi) 

= -^RbMN-BzCl - <JRb)Mi. (8) 

The nuclear Zeeman term is diagonal because the g- 
factors of the two nuclei are equal and not because of 
nuclear exchange symmetry. The N = block of the 
molecular Hamiltonian for a heteronuclear dimer with 
two identical nuclear g-factors would also be diagonal. 

The conservation of the total nuclear spin / and non- 
conservation of M11 and M/2 at high fields may have 
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FIG. 1: Zeeman splitting of hyperfine levels for N = states 
of 85 Rb2 (upper panel) and 87 Rt)2 (lower panel). 



important consequences for the selection rules in spec- 
troscopic transitions used to produce ultracold molecules 
and for the collisional stability of molecules in excited 
hyperfine states. 



B. Zeeman splitting for N — 1 homonuclear alkali 
dimers 



Ultracold homonuclear molecules in N — 1 states arc 
particularly interesting because they are likely to be sta- 
ble with respect to inelastic collisions to produce N = 0, 
at least for collisions with non-magnetic species such 
as other molecules in 1 S states. Such collisions can- 
not change the nuclear spin symmetry and thus can- 
not change N from odd to even. Inelastic collisions 
may well be stronger for collisions of molecules in triplet 
states, because of magnetic interactions between electron 
and nuclear spins. Transitions between odd and even 
rotational levels are permitted in atom-exchange colli- 
sions, such as occur in collisions with alkali metal atoms 

Homonuclear molecules do not possess electric dipole 
moments but do have quadrupole moments. The 



quadrupole-quadrupole interaction is anisotropic and is 
proportional to i?~ 5 , so is longer-range than the R~ 6 dis- 
persion interaction that acts between neutral atoms and 
molecules. The quadrupole-quadrupole interaction aver- 
ages to zero for rotationless states (N = 0), but not for 
N > 0. Quantum gases of rotating homonuclear mole- 
cules may thus exhibit anisotropic effects. 

For N > 0, all the terms in the Hamiltonian (Q]) have 
matrix elements diagonal in N. Some of these are nondi- 
agonal in hyperfine quantum numbers, so the energy level 
patterns are much more intricate. The zero-field splitting 
is dominated in most cases by the electric quadrupole 
interaction and the scalar nuclear spin-spin term. The 
remaining constants (c\ and C3) make much smaller con- 
tributions except for the two L12 isotopomers. For 6 Li2, 
all the terms in the hyperfine hamiltonian contribute sig- 
nificantly. For 7 Li2, the splitting is dominated by the 
electric quadrupole interaction but contributions from all 
the remaining terms are significant. 

Figure [2] shows the "building-up" of the zero-field 
N = 1 hyperfine energy levels for 85 Rb2 and 87 Rb2 in 
three steps: first, only the rotational and the electric 
quadrupole terms are considered; secondly, the scalar 
spin-spin interaction is included; and thirdly, the spin- 
rotation and the tensor spin-spin interaction terms are 
added to complete the hyperfine Hamiltonian. For 85 Rb2, 
the electric quadrupole term alone determines the energy 
level pattern, while for 87 Rb2 there is a significant addi- 
tional contribution from the scalar spin-spin interaction, 
attributable to the relatively large value of C4 for this 
molecule (see tabic fl}. 

The quantum numbers that label the zero-field energy 
levels are included in figure The total angular mo- 
mentum quantum number F is always a good quantum 
number at zero field. In some cases, when there is only 
one pair of values I and N that can couple to give the re- 
sultant F, I is also a good quantum number. Otherwise, 
/ is mixed and the values given in figure[2]are ordered ac- 
cording to their contribution to the eigenstate: the first 
quantum number listed identifies the largest contribu- 
tion. 

The Zeeman splittings for N = 1 states of 85 Rb2 and 
87 Rb2 for different ranges of magnetic fields are shown in 
figures [3] and 0J Each zero-field level splits into 2F + 1 
states with different projection quantum numbers Mp. 
Although in principle both the nuclear (-Hiz) and the ro- 
tational (i?Nz) Zeeman terms contribute to the splitting, 
<?rs.b 9r so that the rotational Zeeman term contributes 
only about 1% for 85 Rb 2 and less than 0.5% for 87 Rb 2 . 

In contrast with the N = case, the Hamiltonian for 
N = 1 is not diagonal and energy levels corresponding 
to the same Mp value display avoided crossings. The 
magnetic field values at which the avoided crossings are 
found, between and 2000 G for 85 Rb2 (lower panel of 
figured]) an d between and 200 G for 87 Rb2 (lower panel 
of figure [3]), scale with the ratio between the electric 
quadrupole constant and the nuclear g-factor. 

For larger magnetic fields, Mj and Mn become indi- 
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FIG. 2: Zero-field hyperfine splitting for N = 1 states of 
85 Rb2 (upper panel) and 87 Rb2 (lower panel). For each 
species, the hyperfine energy levels obtained when only the ro- 
tational and electric quadrupole terms are included are show 
in the left column. The effect of adding the scalar spin-spin 
interaction is displayed in the central column and, finally, the 
right column shows the splitting when the whole hyperfine 
Hamiltonian is considered. All the energies are referred to 
the lowest hyperfine level for the complete Hamiltonian. 
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FIG. 3: Zeeman splitting for N = 1 states of 85 Rb2 (upper 
panel) and 87 Rb2 (lower panel). The inset for 87 Rb2 shows 
the avoided crossing of Mf = states. 



vidually good quantum numbers and the energy levels 
corresponding to the same value of Mj gather together. 
Both features are illustrated in figure |4] where, for the 
sake of clarity, the values of Mm are included only in the 
lower panel. Equation lA. Ill shows that the matrix repre- 
sentation of the nuclear Zeeman term in the spin-coupled 
basis is diagonal with nonzero elements proportional to 
Mi and independent of any other quantum number. As 
the magnetic field increases the nuclear Zeeman terms 
becomes dominant and the slope of the energy levels is 
determined by Mj. 

The results in figure [4] neglect the diamagnetic Zee- 
man interaction, which is not completely negligible at 
the highest fields considered (up to 4000 G). The justifi- 
cation for this is as follows. The Hamiltonian for the dia- 
magnetic Zeeman interaction [l2| consists of two terms 
proportional to the square of the magnetic field: one de- 
pending on the trace of the magnetizability tensor and 
the other is proportional to its anisotropic part. The first 



term has a value around 200 kHz at 4000 G for 85 Rb 2 . 
Although this quantity is not negligible, it has not been 
included because it simply shifts all the energy levels by 
the same amount and has no effect on splittings. The 
second term is diagonal in the spin-coupled basis set and 
its nonzero elements depend on N and Mm- For 85 Rb2 
at 4000 G it would shift the energy levels by about 15 
kHz. It is therefore very small compared to the nuclear 
Zeeman effect. 



V. CONCLUSIONS 

We have explored the hyperfine energy levels and Zee- 
man splitting patterns for low-lying rovibrational states 
of homonuclear alkali-metal dimers in : S states. We have 
calculated the nuclear hyperfine coupling constants for all 
common isotopic species of the homonuclear dimers from 
Li 2 to Cs 2 and explored the energy level patterns in detail 
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90 Hz for 41 K 2 and 160 kHz for 133 Cs2- When a magnetic 
field is applied, each level splits into 21+1 components 
but all the levels with a particular value of Mj are paral- 
lel. This is different from the heteronuclear case, and for 
homonuclear molecules I remains a good quantum num- 
ber in a magnetic field. However, the projection quantum 
numbers M/i and M/ 2 for individual nuclei do not be- 
come nearly good quantum numbers at high fields for 
homonuclear molecules. These differences in quantum 
numbers may have important consequences for spectro- 
scopic selection rules and for the collisional stability of 
molecules in excited hyperfine states. 

Molecules in excited rotational states are also of con- 
siderable interest. In particular, molecules in N = 1 
states may be collisionally stable because transitions be- 
tween even and odd rotational levels require a change 
in nuclear exchange symmetry. Molecules in excited ro- 
tational states have anisotropic quadrupole-quadrupolc 
interactions that are longer-range than dispersion inter- 
actions. The hyperfine energy level patterns are consid- 
erably more complicated for N = 1 states than for TV = 
states and / is not in general a good quantum number 
even at zero field. 

The results of the present paper will be important 
in studies that produce ultracold molecules in low-lying 
rovibrational levels, where it is important to understand 
and control the population of molecules in different hy- 
perfine states. 

FIG. 4: Zeeman splitting for N = 1 states of S5 Rb 2 . The 
avoided crossings for Mf = are shown in the lower panel. 
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APPENDIX: A 



Explicit expressions for the matrix elements of the molecular Hamiltonian terms are now provided. The equations 
are valid for homonuclear molecules. 

The matrix elements for the rotational term (H rot ) are given by 



(™ jv (/i/i)JM J |fl IDt |JV 'M N ,{hh)I> 'M r ) = S NN , 5 Mn m n , &ii> 5 MiMi ,BN(N + 1) (A.l) 
{N{hh)IFM F \H sc \N'{hh)I'F'M F ,) = 5 NN , S ir 5 FF , 5 Mf m f ,BN(N + 1). (A.2) 




The matrix elements for the electric quadrupole interaction (77q) are given by 
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(NMn (7i 7i )IMj I 77 Q \N' M N ,{hh)l' M r ) = (-l) 2/l {4(2iV + 1)(2JV' + 1)(2/+ 1)(2/' + l)} 1 / 2 



AT 2 AT' 




Ji 2 h 
-h /j 



x E 

FM F 



(-1) 



F+2M F 



(2F+1) 



N I F 
M N Mi —Mp 



n' r f 

M N < M r -M F 



INF 

n' r 2 



(A.3) 



(N(hh)IFM F \ H Q \N'{hh)I'F'M F ,) =6 ff ,5 Mf m f , 

' N 2 N' 



(eqQ) 



i(-l) 



JV+JV'+F+2Ji 



x {(2/+l)(2/' + l)} 1 / 2 







Ji 2 /! 
-/i h 



{4(2AT + l)(2Af' + l)} 1/2 



{' 


N 
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5} 


\ N' 


r 






h 





(A.4) 



The matrix elements for the spin-rotation interaction (i?u) are given by 



(NMnihhVMjl H u \N'M N ,(I 1 h)I , M I ,) = 6 NN ,6 II ,c 1 (-l) N - M »- M <+ 2h + 1 2(2I+l) 
x {(2iV + l)N(N + 1)(27 1 + !)/!(/! + l)} 1 ^ J ^ 1 ^ 1 



E 



(-If 



AT 1 AT 
-Mjv p M N > 



I 1 7 

-Mj -p M F 



(A.5) 



{N(hh)IFM F \ 77„ iN'ihl^I'F'Mp,) ^ -26 NN , 5 IF 6 FF , S MF M Fl c 1 (-l) F+2h - MF - I (2F + 1) 
x (2/ + !){(27V+l)iV(7V+l)(2/ 1 + !)/!(/! + l)} 1 /2 1 ^ ' ^ 1 



x E 

F",M F „ L 



E 



(-1) F + m p"(2F" + 1) 

(-iff f 1 f " V 

v y I -M F p M F » / 



N F I 
F" N 1 
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(A.6) 



The matrix elements for the scalar nuclear spin-spin interaction (77 sc ) are given by 



(ArM Ar (7 1 7 1 )7M 7 |77 sc |Ar'M JV ,(7 1 7 1 )7 , M 7 ,} = 5 NN ,8 MN M Nl 5 ir 8 MlMi , -c 4 [7(7 + 1) - 2h(h + 1)] (A.7) 



{N{hh)IFM F \H tc \N'{hh)I'F'M F ,) 



W 8 FF i 5 Mf m f , ^c 4 [I(I + 1) - 27! (7! + 1)] (A.8) 
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The matrix elements for the tensor nuclear spin-spin interaction (H t ) are given by 

(NMn (7i 7i )IMj I H t \N'M N ,{hh)l'M r ) = -cz^fm{-iy- Ml - MN h{h + l){2h + l){{2N +1){2N' + 

1/2 ( N 2 N> 



1)}V» 



x {(2/+l)(2/' + l)} 1 







/! /! 1 
/l /l 1 
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Af 2 JV' W / 2 1' 
-M N p M N , J y —Mi -p M r / 

{N(hh)IFM F \ H t \N'(hh)I'F'M F ,) = -S FF , 5 Mf m fI c 3 V30{-1) n ' +n+i+f h(h + l)(2Ji + 1) 

' h h 1 
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The matrix elements for the nuclear Zeeman term (i?iz) are given by 

(iVMjv(/i/i)/M 7 | i?iz iN'Mwihl^I'Mp) = -5 NN ^M N M N ,SirS Ml M Il gi^B z (l - (Ti)M/ 

(N(hh)IFM F \ H 1Z \N'{hh)I'F'M F ,) - <5 Mf m f ,5iMn5 z (1 - ^i)(-1) 2F " Mf 

x (-l)^ 2 '! (21 + 1){4(2F' + 1)(2F + l)(2/i + + 1)} 1/2 
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The matrix elements for the rotational Zeeman effect (i?Nz) are given by 

{NM N {hh)IM!\ H NZ \N l M N ,(I 1 I 1 )I l M I ,) = -5 NN ,8 MN M N ,5ir5M I M I ,g,mBzM N 
{N(hh)IFM F \ H NZ \N'(hh)l'F'M F ,) = -S NN , S ir 5 MF M F ,g rf i N B z (-l) N+I + F + F '- MF+1 

x {(2F' + l)(2F+l)(2Ar + l)iV(iV+l)} 1/2 
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